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Abstract 

In dimension 1, we show that the Taylor expansion of a "generic" poten- 
tial near a non degenerate critical point can be recovered from the knowl- 
edge of the semi-classical spectrum of the associated Schrodinger operator 
near the corresponding critical value. Contrary to the work of previous 
authors, we do not assume that the potential is even. The classical 
Birkhoff normal form does not contain enough information to determine 
the potential, but the quantum Birkhoff normal form doefl 

In a companion paper |5j, the first author shows how the potential itself 
is, without any analyticity assumption and under some mild genericity 
hypotheses, determined by the semi-classical spectrum. 

1 Introduction 



In this papecl, we will only consider a configuration space of dimension 1. 
Let us consider a (classical) Hamiltonian 
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with V(0) = E , V'(0) = 0, V"(0) = ±1(0). We have 

oo 
3=3 

with n ± = 1(£ 2 ± x 2 ). The Hamiltonian H can be quantized Schrddine;er 
operator H = — |^ 2 J^- + V{x) where the Taylor expansion of V at x = is 
Eq + Ylj^=2 a j x '' w ith a2 = ±2- This operator admits a semi-classical Birkhoff 
normal form [8] (denoted the QBNF) at the origin of which the Weyl symbol is 
a formal power series of the form 

b = n± + b jA 23Qk ± ■ (!) 

2j+k>2 

In this paper, we are interested in the following "inverse spectral problem" : 

does the QBNF, given in (TjQ), of the Schrodinger operator determine 
the Taylor series of V? 

We cannot hope for a positive answer , because V(x) and V(—x) give the 
same QBNF. Moreover 

Remark 1.1 the classical BNF does not suffice to determine the Taylor 
expansion of V at x=0 

Let y = f(x) = x + 0(x 2 ) be an analytic function whose local inverse near 
is of the form x — y + g(y) with g an even function. Then the Hamiltonian 
Hf = 2 (£ 2 + f(x) 2 ) is classically conjugate to Q + near the origin, in particular 
all its trajectories are of period 2ir: it is enough to show that the action integrals 
1(E) = f^2 + ft x \2 <2 E^ x ^ are ^ e same 'i using the change of variable x = y + g(y), 
we get 1(E) = j^ 2+y2<2E (l + g' (y))dydC, and using the fact that g' is odd we get the 

result. A simple example is V(x) = | (\/l + 2x — l) 2 . This result is reminiscent 
of the well known result for Zoll surfaces in Riemannian geometry JMj. 

However, an even potential can be determined by the classical BNF, as a 
consequence of a result of N. Abel F7]@. 

Our main result is: 

Theorem 1.1 The coefficients ±a 3 and are determined from 60,2 one? fr^o by 
the formulas: 

2 5 
a-i = ±V&1,0) «4 = ~bo,2 + -bi, • 

If 03 does not vanish, all dj 's are determined from the 6o,fe 's and the b^k 's 
once we have chosen the sign of a%. 

3 Assuming that V"(0) — ±1 does not affect the results below, because 02 = V"(0)/2 is 
known from the first eigenvalue if 02 > and from the density of states if 02 < 
4 We are grateful to Hans for pointing this out to us 
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This result is reminiscent of the much more sophisticated results by Zelditch on 
the Kac problem [9j. If we use a (trivial) particular case of the result of [6], we 
get 

Corollary 1.1 If we know the asymptotic expansions of the eigenvalues X n (h) for 
all n's of a Schrddinger operator near the minimum x = of the potential and 
V"(0) > 0, we know the value ofV'"(0) and, ifV"'(0) ^ 0, the Taylor expansion 
of the potential at that point. 

In fact, we have the more precise result: 

Corollary 1.2 From the knowledge of the N first eigenvalues of H modulo 0(h 2N ), 
one can recover the Taylor expansion of V to order 2N. 

A similar result holds for a local non degenerate maximum of V using the 
"density of states" techniques. This is the content of Section [TUl 

Corollary 1.3 If E is a non degenerate local maximum of V and is the only 
critical point ofV on the set V = E , the knowledge of the semi- classical spectrum 
of H in some intervall }E , Ei[ (or ]E 1: E [) determines V"'(0) and, provided that 
V'"(0) 7^ 0, the Taylor expansion ofV at x = 0. 

and also in the case of a local minimum (Section I10.4p : 

Corollary 1.4 If E is a non degenerate local minimum ofV and is the only 
critical point ofV on the set V = Eq, the knowledge of the semi- classical spectrum 
of H in some interval }Ei,E2[, with E\ < Eq < E%, determines V"'(0) and, 
provided that V"'(0) ^ 0, the Taylor expansion ofV at x = 0. 

Knowing the semi-classical spectrum as a function of h~ seems to be an huge 
amount of information. As was showed in [3] , this is however the case for the ef- 
fective Hamiltonians driving the propagation of waves inside a stratified medium. 

2 A counterexample for a general Hamiltonian 

The QBNF of a general Hamiltonian, independent of h, H(x,C,) = Q± + 0(3) 
is not enough to know the Taylor expansion of H at the singular point. It is 
enough to consider H — | ((£ — 3x 2 ) 2 + x 2 ) which is gauge equivalent to Vt + by 
the gauge transform u — > ue lx3 . 

3 Review of the Moyal product 

The Moyal product is the product rule of symbols of Weyl quantized \1/Z)0's, it 
is given by: 
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with 



{a, b} 3 == E ( J ) (-l) P ^r Pa ^" P ^ • 



p=0 

We will also use the Moyal bracket, 

[a,b]* — a*b — b* a . 

We have 



i ~ i f h\ 2j 

In particular, {a, b}i = a^b x — a x b^ is the Poisson bracket and 
{a, b} 3 = a^b xxx — 3a^ x b xx ^ + 3a 5 

We have: 

% h 2 h 4 

-[a, b]* = {a, 6}i - —{a, b} 3 + j^{a, b} 5 + ■ ■ ■ . 

4 The Weyl algebra 

The 'Weyl algebra" which consists of formal power series in h and (x, £ ) 

oo 

where W^- is the space of polynomials in (x, £ ) and h of total degree j and the 
degree of x l ^ m h n is / + m + 2ra. 14 7 is a graded algebra for the Moyal product: we 
have Wj * W k C H^+fc and hence W^]* C W j+k -2- Moreover, if we define 

Wj~ as the subspace of Wj which is generated by monomials of even degree in h, 
we have 

l j.[w;,w k + ]*cw; +k _ 2 ; 

we will define W + = Wj~ which is a Lie algebra for the bracket | [., .]*. W + 

is the (formal) Lie algebra of FIO's which are tangent to the identity at the the 
origin. The grading is obtained by looking at the action on the (graded) vector 
space of symplectic spinors: if F = YljLo ^ Fj{X) with F e «S(R), we define 
fn{ x ) — h~^F{x/K) whose microsupport is the origin. W + acts on this space of 
functions in a graded way as differential operators of infinite degree: if w e W, 
w.f = OP h (w)(f). 
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5 Moyal versus functional QBNF 



There are two different QBNF: 

• The first one is a Weyl symbol B = ^2 bj :k h 2j £l k as before, 

• The second one is an operator which is a formal power series of the harmonic 
oscillator Cl of the form B = ^ bj ! j e h 2 ^Cl k . 

The second one is the Weyl quantization of the first. So they are equivalent. 
The equivalence can be made explicit in both direction by computing Op Weyl (f2 fc ) 
or the Weyl symbol of tl k . The functional form is useful in order to compute 
successive approximations of the eigenvalues, while the Weyl form is easier to 
compute using the Moyal product. 

6 Useful Lemmas 

The following result is classical: 

Lemma 6.1 The equation {Q±,P}i = Q where Q is a given homogeneous poly- 
nomial of degree N has a solution P , a homogeneous polynomial of degree N , 

• if N is odd 

• if N = 2N' is even and c±(Q) = where c± is a linear form on the space 
of homogeneous polynomials of degree N which satisfies c±(Q± ) = 1. In 
particular, given Q, the equation {Q±,P}i = Q — c±(Q)Q± has a solution. 

Remark 6.1 In the case c + (Q)£l+ is the average of Q under the natural 
action of S 1 on homogeneous polynomials of degree 2N' . 

Definition 6.1 We will denote by S^ v _ 1 the homogeneous polynomial of degree 
2N — 1 which satisfies 

{fl±, S 2A r_ 1 }i = X 

Lemma 6.2 We have 
We can also check the: 

Lemma 6.3 The polynomials x 2N ' are not Poisson brackets of the form x 2N ' = 
{Ot,P}i, i.e. c ± (x 2N ')^0. 
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7 The QBNF 



In order to reduce to the QBNF, we will use automorphisms of W + of the form 

H ^H s = exp(iS/h) * H * exp(-iS/h) = exp ^ad(S)*^ if 
with S* = S*3 + 5*4 + • • • G W+. We get: 

^ k brackets 

ff s = ff + jM* + -. + i(lj '[5, [^••,[s,#]T---]* + --- , 

which is a convergent formal power series whose k— th term is of degree > fc + 2. 
The brackets will be calculated using the Moyal bracket. We remark that the 
terms of degree in h give the calculation of the classical BNF (denoted CBNF) 
where the brackets are now just Poisson brackets. 



8 The first terms 

Let us consider V(x) = \x 2 + ax 3 + bx A -\ whose QBNF is Q+An 2 +BH 2 +0(6) 

where 0(6) means terms of degree > 6 in the Weyl algebra. Here we assume 
= \(i 2 + x 2 ). Our first result is: 

Theorem 8.1 

A = — -a 2 + -b, B — a . 
4 2 

The calculation: we start with S = S 3 + S4 where S 3 (x, £)(resp. S4 (#,£)) is a 
homogeneous polynomial of degree 3(resp. 4). There is no need to put terms in 
h 2 in £4 because they would be of the form cti 2 which is in the center. We have 
then: 

exp(i[S, ]*)H = n+ l -[S,H}*+ 1 - (£) 2 [S, [S, H]r + 0(6) . 

By identification of terms of degree 3 and 4 and using the expression of the Moyal 
bracket [.,.]*: 

l\f,9]* = {f,9h-^t?{f,9h + -- - , 
we get the system of equations: 

j (3) ax 3 + {S 3 ,Q} 1 =0 

\ (4) bx' l + {S 3 ,ax 3 } 1 + {S i ,n} 1 + l{S 3 ,{S 3 ,n} 1 } 1 -j- 4 h 2 {S 3 ,ax 3 } 3 =An 2 + Bh 2 

Using Equation (3), Equation (4) splits into 2 equations: 

/ (4') -±{S 3 ,ax 3 } 3 =B 

\ (4") bx A + \{S 3 , ax 3 }! + {S 4 , = AVI 2 
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From Equation (3) and the formula for £2iv-i given in Lemma [6.11 we get 

S 3 = -a(x 2 e + ^ 3 ) • (2) 
From Equation (4'), we get B = a 2 . From Equation (4"), we get the value of A. 

9 The induction 

We carry out the proof in the case ofQ + and Eq — 0. The minus case is similar. 
Let us start with 

H' = tt + + a 3 x 3 + ■■■ + a 2N -2X 2N ~ 2 
and S' = S 3 + S4 + ■ • • + S2N-2 with Sj G Wj, so that 

exp U[S', A H' = Q + + B A + ■ • • + £ 2iV _ 2 + R 2N _, + R 2N + ■ ■ • (:= B') , 

with 

• B 2 j G W 2 + j a polynomial in h 2 and Q + 

• For n = 2N - 1 and n = 2N, R n G W n . 

In other words 5" generates the transformation which converts H' into its QBNF 
mod 0(2N — 1). The polynomials H' and S' and the partial QBNF B' are known 
by the induction hypothesis. We are now trying to get S" = S 2 n-i + S 2 n so that 
S = S' + S" converts H = H' + ax 2N ~ l + bx 2N into the QBNF mod 0(2N + 1). 
We will only consider the terms of degree and 2 in H. So we can split every 
polynomial Pj in into Pj = P® + h 2 P 2 + ■ ■ ■ with P 2 of degree j — 4 in (x, £). 
The equation to solve is: 

exp (US' + S", .]*) {H' + ax 2N - 1 + bx 2N ) = 

n + + B 4 + --- + B 2N „ 2 + B 2N + 0(2N + 1) [ } 

with B 2N = &2v^+ + &2V^- 2 ^+ -2 + ' ' ' • We hope to recover a and b from 6°^ and 
b\ N using what we know already at this step. The left handside of Equation (j3J) 
splits into: 

exp (US' + S", A (ax^+bx™) = ax^-'+bx 2 " + % -{S 3 , ax 2N - 1 Y+0(2N+l) 
and 

exp U[S' + S", Aj H' = B'+ % -{S", tt + + a 3 x 3 }* + 

1 / • \ 2 
III 



+ 2 U ) (t^iv-i, [5 3 , n+fi* + «+]*]*) + 0(2N + 1) . 



So that, we get 
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• In degree 2N — 1: 

ax ™-l + {S O n _^ + R N i = 

{S 2 jy_l, + R-2N-1 = • 

We see that S , f JV _ 1 is known at this step, while S® N _ 1 is modulo known 
terms a solution of 

{fi + ,s 2 Vi}i = ^ 2A ^ • 

This equation gives, always mod known terms: 

with S 2 jv-i given by Definition 16.11 

• In degree 27V: 

bx 2N + l - ([S 3 , ax 2N - 1 }* + [S 2N , n + f + [S 2N ^, a 3 x*Y) + 

+\ (0 ([S2N-1, [s 3 , + [s 3 , [sW-i, + = b 2N + o(2iv + 1) . 

The previous equation gives one equation in h° and one in h 2 : 

— degree 2N, hbar 

bx 2N + {S 3 , ax 2 ^ 1 }! + {S° 2N , + {S 2 Vi, a 3 x 3 }! + 

+ 2{^27V-l) {£3, ^+}l}l + 2 i^ 3 ' {&2N-li ^+}l}l + ^2iV = &2AT^+ 

This can be simplified as: 

{^+,S 2Ar }i = ~b 2 N^++bx 2N +-{S 3 ,x 2N 1 }i+-^-{S 2N _ 1 ,x 3 }i+R2 N —-{S 3 ,R2 N _ 1 }i 

(4) 

This gives fe 2Ar = ^6+7^030 modulo known terms. Moreover, Lemma 
16.31 implies (3n 7^ 0. 

— degree 2N ', hbar 2 

-1 ({S 3 , ax^^ja + a 3 x%) + {Sf„, + 

+ 2 ({^2iV-l> {S3, ^+}l}l + {S3, {S 2 ^.!, fi+}i}i) - 
~48 ("t^JV-D {S3, ^+}l}3 + {S3, {S2JV-I, ^+}l}3) + Rm = &2iV^+ 2 



S 



which can be simplified as: 

{Q + , S% N }i = —bl N Q+ 2 ——[a{S 3 ,x 2N 1 } 3 + a 3 {S% N -i, x 3 } 3 )+Rl N , 

^ (5) 

modulo known terms. This gives, using Lemma I6.1[ b\ N = 5nciQ>3 
modulo known terms. 

From Equation (jSJ) and the expressions for 5*3 (Equation (j5J)) and T,2N-i 
(Lemma l6.2p . we get: 

m c2 i (N-l)(2N 2 -4N + 3) 2N _ A , 2 ^ N _ 2 

o 2 n\i — aa 3 g £ — "2Af"+ mod known terms 

(6) 

Because x 2N 4 is not a Poisson bracket with Q + by Lemma |6.3[ we get 6~n ^ 0. 

From the fact that (3^ and 5n do not vanish, this concludes the induction 
N-l^N. 



10 Getting the QBNF from the density of states 
in case of a local extremum of the potential 

10.1 /i- dependent distributions 

Let T h be an ^-dependent Schwartz distribution on an open interval J. 
Definition 10.1 The family is 

• regular at the point E G J if there exists a sequence of functions Tj which 
are smooth in some neighbourhood K of E with j = j , jo + 1, ■ ■ ■ (jo G Zj, 
so that, for any f G C™(K), we have the asymptotic expansion T n (f) = 
Et= j0 V fjT^fWdx. 

• right regular (resp. left regular ) at the point Eq G J if there exists E\ > Eq 
(resp. Ei < Eq) and a sequence of functions Tj which are smooth in 
some neighbourhood of Eq with j = j , jo + 1, • • • (jo G Z,), so that, 
for any f G C%°(]Eo, Ei[), we have the asymptotic expansion T^(/) = 

SMx)f{x)dx. 

We will use the following notations: 

Definition 10.2 If T^ is a family of distributions on J and E G J, T£ (resp. 
T h ), the right (resp left) singular part of T n is the equivalence class ofT h modulo 
families of distributions which are right- (resp. left-)regular at the point E . 
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10.2 Density of states 

Consider a smooth potential V : / — > R where 7 is an open interval with & I 
and liminf x ^a/ = > — oo and let H be the Schrodinger operator with 
potential V. 

Definition 10.3 The density of states is the h— dependent Schwartz distribution 
Tn on] — oo, Eoo[ defined by 

D h (f) := Trace/ (tf ) . 

Lemma 10.1 If J is an open subset of] — oo, -Eoo[ which contains no critical 
values ofV, the density of states is regular at every point of J. 

Proof. - 

Let us denote by H = |£ 2 + V(x) the symbol of the Schrodinger op- 
erator H. The operator f(H) is a pseudo-differential operator whose 
symbol f*(H) is given (see [1]) by: 

j>i, i>i 

where the Pj/s are smooth functions locally computable from the 
symbol H. It is now enough to check that / — > (27r^) _1 J J Pjj(x, £) f( l \H)dxd£ 
is regular at each point of J using the fact that H has no critical value 
in J. 

□ 

10.3 Singularity of the density of states near a local max- 
imum of the potential 

Let us assume that V(0) = E < E^, V'(0) = and V"(0) < 0. Assume also 
that is the unique critical point of V whose critical value is E Q . 
We have the: 

Theorem 10.1 If the QBNF of H is 

2j+k>2 

the density of states is right and left singular at the point E Q and one can recover 
the full QBNF (the coefficients bj^) from the right (resp. left) singular part 
(resp. D^) of the density of states D h at E . 

In what follows, it is more convenient to use f2_ = x£. 
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10.3.1 The singularity of the density of states and the QBNF 

Lemma 10.2 If B is the QBNF of H, the singular part of the density of states 
is the same as that of the family of distributions 

where f*(B) is the Weyl symbol of f(B) and D is the square max(|x|, |£|) < 1. 
Proof. - 

Let II = Op Weyl (cj) be a compactly supported *$fDO whose Weyl 
symbol uo is = 1 near (0, 0). We have 

D h (f) = (27rh)- 1 (J Jou*f*(H)dxd£+ J J \l-u) *f*(H)dxd£ 

Using (the proof of) Lemma 110.11 the second term is a regular dis- 
tribution. The first term can be transformed using the QBNF: there 
exists an FIO U, microlocally unitary, which transforms H into its 
QBNF and hence for every function /, we have: 

V*f(B)U = f(H) 

microlocally near the origin. In this way, we get 

Trace(n o f(H)) = Trace (IK/*/ (B)U) . 

Introducing IL := UHU* (a ^>DO whose Weyl symbol is = 1 near 
the origin) and using the commutativity of the trace, we have: 

Trace (n of(H)) = Trace (IT o /(£))) . 

It remains to check that, if IT = Op Weyl (cui), / — > J f, IxR \\ D 0Ji * 
f*(B)dxd£, is regular. 



10.3.2 Computing some singularities 

Lemma 10.3 Let us consider the family of distributions 



□ 
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on }E ,Ei[ (we consider only the case E\ > E , the other case is similar), where 
the bj 's are smooth on [E , Ei] and b (u) = E Q + Y^jLiflj u ^ A > 0. Then 
Kh{f) admits an asymptotic expansion in powers of h: 

oo 

K h {f) = Y. K M) h2j 

j=0 

and the right singularities of K , - • ■ at the point E determine the Taylor 

expansions of bo, - ■ ■ , b^ at the origin. 

Proof. - 

Let us Taylor expand the integrand as: 

/ (£r=o ^M)) = f(bo«)) + f'(boM) (Er=i MM)) 

= f(b (xO) +f> 2j (fMx^bM) +j2f {l) ( b MmA^ , 

where the functions Rjj depend only on 61, • • • , bj-i. 
We have to prove the following 2 facts: 

I. The right singularity of 



J Jj{h{xi))dxdi 



determines the Taylor expansion of 60 at the origin. 
2. The right singularity of 



J J D f'(b (xO)bM)dxdC 



determines the Taylor expansion of bj at the origin, assuming 
the Taylor expansion of b is known. 

Both are easy consequences of the following elementary calculus re- 
sult: 



/ / f(bo(xO)bM)dxd^= [ V(%(co(0K(0|l°g(*-£o)|^ 
J Jd Je 

(modulo smooth distributions) where cq is the inverse function of bo- 



□ 
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10.3.3 End of the proof of Theorem [TOTTI 

We have 

oo 

nm*>) = E 27T (/ (2j) (5(^o))(5 - sczb))*^) («b) • 

j=0 3 ' 

It is enough to check the: 
Lemma 10.4 If 

oo 

P(B) = f(B) + E ^' E f W ( B )Ru > 
i=i / 

i/ie 's depend only on bo, - ■ ■ , fej-i- 
Proof. - 

The *— powers oi B — B(zq) evaluated at zo start with terms in ft 2 
and the 6;'s, for Z > j have already an ft 2 - 7 in front of them! 

□ 

So everything works as if f*(B) = f(B) and we are reduced to Lemma H 0.31 
10.4 The case of a local minimum 

The same strategy applies, but now the density of states is right AND left regular, 
with a jump singularity at the point E . 
We get: 

Theorem 10.2 If the QBNF of H is 

B = Q + + E b j,kh 2j tt k + , 

2j+k>2 

the density of states is singular at the point E and one can recover the full QBNF 
(the coefficients bj t k) from the singular part of the density of states at E . 

The proof is very similar to the case of a local maximum. We have now 
a "Heaviside singularity" , meaning that the density of states is right AND left 
regular, but the functions Tj defined by 

oo „ 

DM = E / Wit 

i=-i J 

and their derivatives have jumps at the point E Q . We have only to look at the 
singularities of T : / — > J f(Q + )dxd£. We have T(f) = 2ir J °° f(u)du, so 
T = 2ttY where Y is the Heaviside function. 
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11 Open problems 



• Is the result still true if 03 = 0? This is the case modulo some global 
assumption on V (see [5]). In fact in [5], it is shown that, modulo some 
genericity assumptions, the potential itself is determined from its semi- 
classical spectrum. 

• Is the result still valid in any dimension? We think that the answer is 
no, at least it does not work with the same arguments; let us assume that 
the quadratic part of the Hamiltonian is H2 = uiQi + U2Q2 with Qi (resp. 
f^) harmonic oscillators in x\ (resp. 22). 

— Non resonant case: u\ and U2 are independent over Z. In degree 4, 
the QBNF has 4 unknown coefficients, an homogeneous polynomial 
of degree 2 in (f^,^) and the coefficient of h 2 . On the other hand, 
^3(21, a; 2 ) + ^4(21,22) has 9(> 4) coefficients. However, it is possible 
that higher terms in the QBNF give other information's... 

— Resonant case: iO\ = u>2- In degree 4, the classical BNF has already 
9 coefficients (it is a polynomial of degree 4 on M 4 invariant by the 
circle action generated by the flow of Q\ + f^), this seems promising. 
However, we have to take into account an 0(2) action by isometries 
in IR 2 : on one hand, we can only expect to determine the potential up 
to this action; on the other hand, the QBNF is determined only up to 
action by SU(2). 



12 Homogeneity properties of the QBNF 

We have the following: 

Theorem 12.1 The bj^'s (coefficients ofh 2 ^Vt k in the QBNF) satisfy the follow- 
ing homogeneity properties: 



bj,k(ta>3, t 2 a4, ■ • ■ , t n a n+ 2 

Proof. - 

Let us consider 
1 



)=t 2 ^-\ k (a 3 ,a,,...) 



d 2 

-ti 2 -^ + ) + ta 3 x 3 + ■■■ + t n ~ 2 a n x n + ■ 



and make the change of variable tx = y, Hi — t 2 h. We get a new 
operator 

d 2 
dy 2 



t" 



- h liri + y 2 ) + a 3V 3 + ■■■ + a n y n + ■ 
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The spectrum of the second one is then t 2 times that of the first one. 
This implies the property. 



□ 



References 

[1] N. Abel. Auflosung einer mechanichen Aufgabe, Journal de Crelle 1:153-157 
(1826). 

Arthur L. Besse. Manifolds all of whose Geodesies are closed. Springer. 
Ergebnisse no 93 (1978). 

Yves Colin de Verdiere. Mathematical models for passive imaging II: Effec- 
tive Hamiltonians associated to surface waves. \ArXiv:math-ph/0 6 10044- 

Yves Colin de Verdiere. Bohr-Sommerfeld rules to all orders. Ann. Henri 
Pomcare 6:925-936 (2005). 

Yves Colin de Verdiere. A semi-classical inverse problem II: smooth poten- 
tials. Preprint (December 2007). 

Victor Guillemin, Thierry Paul & Alexandro Uribe. "Bottom of the well" 
semi-classical wave trace invariants. ArXiv:math-sp/060861 7 and Math. Res. 
Lett. 14:711-719 (2007). 

Victor Guillemin & Alexandro Uribe. Some inverse spectral results for semi- 
classical Schrodinger operators. ArXiv:math-sp/ 0509290 and Math. Res. 
Lett. 14:623-632 (2007). 

San Vu Ngoc & Laurent Charles. Spectral asymptotics via the Birkhoff 
normal form. \ArXiv:math-s p / '0605096. 

Steve Zelditch. The inverse spectral problem. Surveys in Differential Geom- 
etry IX, 401-467 (2004). 



15 



